Technical  Report  1404 

December  1990 


Detection  of  Targets  in 
Terrain  Clutter  by  Using 
Multispectral  infrared 
Image  Processing 


L.  E.  Hoff 
J.  R.  Evans 
L.  E.  Bunney 


OTIC 

ELECTE  f 
JUL&21S9U 

B 


D 


Approved  for  public  release;  distribution  Is  unlimited. 


NAVAL  OCEAN  SYSTEMS  CENTER 

San  Diego,  California  92152-5000 


J.  D.  FONTANA,  CAPT,  USN  H.  R.  TALKINGTON,  Acting 

Commander  Technical  Director 


ADMINISTRATIVE  INFORMATION 

The  work  reported  here  was  performed  by  members  of  the  Processing  Applications  Branch, 
Space  Systems  and  Technology  Division,  Surveillance  Department,  with  funding  provided  by  the 
Defense  Advanced  Research  Projects  Agency  (Code  AVSTO),  1400  Wilson  Boulevard, 

Arlington,  VA  22209. 


Released  by 

D.  N.  Williams,  Head 

Processing  Applications  Branch 


Under  authority  of 
F.  M.  Tirpak,  Head  (Acting) 
Space  Systems  and 
Technology  Division 


ACKNOWLEDGMENT 

The  authors  would  like  to  thank  Larry  Stotts,  DARPA,  and  Ed  Winter,  Technical  Research 
Associates,  who  have  also  been  deeply  involved  in  the  development  of  signal-processing  algorithms 
for  multispectral  imagery.  This  work  was  partially  funded  by  the  DARPA  Multi-Spectral  Infrared 
Camera  (MUSIC)  program. 


SUMMARY 


PROBLEM 

Investigate  a  weighted-difference  signal-processing  algorithm  for  detecting  ground  targets  by  using 
dual-band  IR  data. 


RESULTS 

Three  variations  of  the  algorithm  were  evaluated:  (1)  simple  difference;  (2)  minimum  noise;  and 
(3)  maximum  SNR.  The  theoretical  performance  was  compared  to  measured  performance  for  two 
scenes  collected  by  the  NASA  TIMS  sensor  over  a  rural  area  near  Adelaide,  Australia,  and  over  a 
wooded  area  near  the  Redstone  Arsenal.  The  theoretical  and  measured  results  agreed  extremely- well. 
For  a  given  correlation  coefficient  and  color  ratio,  the  amount  of  signal-to-noise  ratio  gain  can  be 
predicted.  However,  target  input  SNRs  and  color  ratios  can  vary  considerably.  For  the  targets  and 
scenes  evaluated  here,  the  typical  gains  achieved  ranged  from  a  few  dB  loss  (targets  without  color)  to 
a  maximum  of  approximately  20  dB. 


Accession  For  "T| 


NTIS  GRA4I 
DTIC  TAB 
Unannounced 
Justlf tcatio 

- 

□ 

□ 

n  _  ...  - 

Rv 

Distribution/ 

Availability  Codes 

Dist 

Avail  i 
Spec 

and/or 

lal 

CONTENTS 


1.0  INTRODUCTION .  1 

2.0  MULTISPECTRAL  PROCESSING  CONCEPTS  .  1 

3.0  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS  .  6 

3.1  Experimental  Results  from  the  Adelaide  Scene .  6 

3.2  Experimental  Results  from  the  Redstone  Scene .  14 

4.0  CONCLUSIONS .  18 

REFERENCES .  19 

APPENDICES 

A.  Details  of  Numeric  Results,  Adelaide  Data  Set . A-l 

B.  Model  for  Two-Channel  System . B-l 

C.  Analysis  of  Data,  Adelaide  512  x  512  Scene . C-l 

D.  Details  of  Numeric  Results,  Huntsville  Data  Set . D-l 

FIGURES 

1.  Processing  flow  for  weighted-difference  algorithm .  3 

2.  Weighted-difference  algorithm  variants  for  a  background  band-to-band 

correlation  of  0.9 .  5 

3.  Three  cases  of  relative  target  and  background  temperatures  .  7 

4.  Channel  radiance  versus  temperature .  8 

5.  Adelaide  scene .  9 

6.  Adelaide— weighted-difference  output .  10 

7a-b.  dB  gain  as  a  function  of  color  ratio .  12 

8.  dB  gain  as  a  function  of  input  SNR .  13 

9.  Huntsville,  Alabama— Redstone  scene .  15 

10.  Huntsville,  Alabama— Redstone  2.5-m  tree-line  tanks .  16 

lla-b.  Cumulative  distributions  for  the  unprocessed  input  data  and  preprocessed 

MIN_NOISE  1  and  4  output  data  .  17 

12.  dB  gain  as  a  function  of  color  ratio .  18 

TABLES 

1.  Variants  of  the  weighted-difference  algorithm  .  5 

2.  Detection  statistics  for  52  targets  in  the  Adelaide  scene .  11 

iii 


1.0  INTRODUCTION 


The  detection  of  small  targets  in  clutter  is  a  problem  of  critical  importance  to  wide-area  surveil¬ 
lance  [1].  Stealth  technology  has  made  targets  more  difficult  to  detect.  Targets  can  also  hide  in  clutter 
that  is  either  manmade  or  natural.  One  technology  continuously  being  investigated  for  improving  small- 
target  detection  in  clutter  is  infrared  (IR)  imagery.  Multispectral  IR  films  have  been  effective  for 
detecting  some  forms  of  camouflaged  targets  and  other  targets  using  false  coloring  [2].  Recently,  inter¬ 
est  in  using  multispectral  sensor  data  in  signal  processing  was  stimulated  by  an  analysis  of  such  data 
collected  by  the  NASA  Thermal  Infrared  Multispectral  Scanner  (TIMS)  sensor  in  Australia  during 
1985  [3].  These  data  and  their  analysis  [4]  suggested  that  multispectral  processing  techniques  could 
greatly  improve  detection  of  small,  hard-to-find  targets. 

The  work  presented  here  is  part  of  the  DARPA  Multi-Spectral  Infrared  Camera  (MUSIC)  pro¬ 
gram.  Several  algorithms  have  been  developed  and  tested  on  infrared  data  [5,6,7]. 

The  objective  of  this  work  was  to  investigate  weighted  difference,  a  signal-processing  algorithm  for 
detecting  ground  targets  by  using  dual-band  IR  data.  This  algorithm  is  derived,  analyzed,  and  com¬ 
pared  to  empirically  measured  results. 

Section  2.0  gives  a  derivation  of  the  algorithm  and  its  theoretical  signal-to-noise  ratio  (SNR)  of 
input  to  output.  Section  3.0  gives  the  results  of  applying  the  algorithm  to  two  image  scenes,  and  com¬ 
pares  the  results  to  the  theoretically  predicted  performance  (showing  excellent  agreement).  Numerical 
details  supporting  this  work  have  been  included  in  the  appendices. 

2.0  MULTISPECTRAL  PROCESSING  CONCEPTS 

There  are  two  principles  we  exploit  in  our  use  of  multispectral  IR  imagery.  The  first  is  that  IR 
images  of  natural  backgrounds  in  different  long-wave  infrared  (LWIR)  bands  tend  to  be  highly  corre¬ 
lated.  This  high  degree  of  correlation  allows  the  data  collected  in  one  spectral  band  to  be  used  to 
cancel  the  background  clutter  in  the  data  collected  in  another  spectral  band.  The  second  principle 
involved  is  that  the  LWIR  signatures  of  manmade  objects  generally  do  not  exhibit  the  same  relation¬ 
ship  between  bands  as  do  the  correlated  background  signatures.  We  refer  to  this  difference  in  target 
signature  between  bands  as  “coloring.” 

The  algorithm  for  weighted  differences  is  the  following:  the  scene  data  in  one  band  are  multiplied 
by  a  constant  and  subtracted  from  the  scene  data  in  the  other  spectral  band  on  a  pixel-by-pixel  basis. 
The  constant  is  chosen  so  that  the  background  clutter  in  the  resulting  image  is  minimized.  While  the 
background  is  minimized,  the  target  will  leave  some  residual  in  the  output  image  because  of  its  color¬ 
ing  between  spectral  bands.  Target  coloring  and  background  correlation  determine  the  performance  of 
the  weighted-difference  algorithm. 

The  output  radiance  for  each  band  is  the  sum  of  target  signal  and  background  clutter.  The  back¬ 
ground  clutter  is  assumed  to  be  made  up  primarily  of  natural  vegetation  or  clouds  (which  are  grey 
bodies),  so  the  clutter  will  be  highly  correlated  from  one  band  to  the  next.  Sensor  noise  and 
“colored,”  or  nongrey,  bodies  in  the  image  will  reduce  the  correlation  coefficient  between  the  two 
bands.  The  targets  are  assumed  to  be  manmade  objects  and  are  characterized  as  selective  radiators, 
which  can  have  a  different  radiance  level,  or  color,  in  each  band. 

The  sensor  either  scans  or  uses  a  focal-plane  array  to  collect  data  in  a  rectangular  image.  The 
data  consist  of  pixels  that  are  modeled  as  being  either  target  or  backg  ound  clutter.  Each  pixel 
records  a  measure  of  the  received  IR  radiance  for  that  image  cell.  Th  i  target  pixel  obscures  the  back¬ 
ground,  so  these  pixels  contain  only  target  radiance  and  no  background. 
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Let  the  targerpixel  be  given  by  Si  and  clutter  pixels  be  labeled  as  N-t,  where  i-x  or  y  for  bands 
x  and  y,  respectively.  The  target  radiance  in  the  bands  is  given  by 

Si  =  esiB(T)  i  =  x,y  (1) 

and  the  background  radiance  is  given  by 

Ni  =  eciB(Tc)  i  =  x,y  (2) 


where 


T  -  the  target  temperature 
Tc  =  the  temperature  of  the  background 
eSi,  eCi  =  the  emissivity  coefficients  for  target  and  clutter 
B  =  Planck’s  black-body  function  . 

Ni  is  a  random  variable  because  of  random  variations  from  pixel  to  pixel  of  the  emissivity  and  clutter 
temperature.  Let  mi  and  v;  be  the  mean  and  variance  of  the  clutter  for  band  i.  The  standard  devia¬ 
tion  is  the  square  root  of  the  variance  and  will  be  denoted  as  oj .  The  correlation  coefficient  of  the 
two  bands  is  then  given  by 


Cov(Nx,Ny) 

OxOy 


(3) 


Figure  1  gives  a  block  diagram  for  the  weighted-difference  algorithm.  The  algorithm  multiplies  the 
image  from  band  y  by  a  constant  weight,  ft) ,  and  subtracts  it  from  the  image  from  band  x.  The  objec¬ 
tive  of  the  algorithm  is  to  produce  an  output  image  that  enhances  the  target  signal-to-noise  ratio. 
There  are  three  variants  of  the  algorithm,  depending  on  how  the  constant,  ft) ,  is  derived. 

Variant  1:  Max_SNR 


Choose  ft)  to  maximize  the  output  SNR. 


Variant  2:  Min  Noise 


Choose  ft)  to  minimize  the  variance  of  the  output  image. 


Variant  3:  Sim  Diff 


Choose  ft)  =  1.0. 

Define  the  input  pixel  signal-to-noise  ratio  to  be 

psm®  =  -.m  . 


(4) 
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Figure  1.  Processing  flow  for  weighted-difference  algorithm. 


The  signal-to-noise  ratio  is  measured  in  dB,  given  by  the  following 

SNR/n  (dB)  =  10  log  [psnr(i)] 2  (5) 

where  /  denotes  the  band. 

Let  V  be  the  output  variance  when  only  clutter  is  present;  then  the  output  signal-to-noise  ratio  that 
is  to  be  optimized  is  SNR0U(,  where 


SNR0„, 


(■S’*  ~  (OSy)2 

V 


The  output  variance  can  be  shown  to  be  a  function  of  the  input  statistics  as 


(6) 


V  =  dx  -  lapOxOy  +  ft)2Oy  . 


(7) 


The  three  variants  of  the  weighted-difference  algorithm  can  be  derived  from  the  above  equation. 
However,  the  derivation  is  simpler  and  easier  to  implement  if  each  channel  is  normalized  to  have 
mean  equal  to  zero  and  variance  equal  to  one.  The  output  SNR  becomes 


SNR0«f 


(SnX  ~  (0Sny)2 
(1  -  2cop  +  to2) 


(8) 
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where 


(5/  —  7»;) 
o H 


Sni  — 


.i  =  x,y 


(9) 


The  weiglit  for  variant  1,  <iWx_snr,  can  be  determined  by  taking  the  derivative  of  the  SNR0U*  with 
respect  to  co  and  setting  it  to  zero  to  find  the  critical  points.  There  are  two  critical  values  (see  appen¬ 
dix  B).  The  one  that  maximizes  the  SNR0Uf  is  variant  1: 


^r.iax_SNR 


P  “  c 
1  -  pc 


where  c  is  defined  as  the  ratio  of  the  pixel  signal-to-noise  ratios. 


(10) 


psnrm-m 
psnr max 


-  1  <  C  <  1  . 


(11) 


When  c  =  p  the  weight  is  zero.  This  can  be  interpreted  to  mean  that,  for  the  Max_SNR  algorithm, 
the  two  bands  should  not  be  combined  and  the  best  output  image  when  p  =  c  is  the  input  image  from 
band  x.  Note  also  that,  when  c  =  1/p,  the  weight  is  infinity.  This  can  similarly  be  interpreted  to  mean 
that,  for  the  Max_SNR  algorithm,  the  two  bands  should  not  be  combined  and  the  best  output  image 
when  c  =  1/p  is  the  input  image  from  band  y.  To  avoid  potentially  large  values  of  co  ,  the  input  bands 
were  arranged  so  that  band  x  was  always  the  band  with  the  largest  psnr. 

To  implement  the  Max_SNR  algorithm,  we  need  to  know  the  color  ratio  of  the  target  pixels.  For 
selective  radiators,  this  ratio  would  not  be  known  ahead  of  time.  For  the  calculations  here,  we  knew 
where  the  target  pixels  were  and  measured  the  color  ratio  from  the  actual  data,  which  is  not  a  very 
practical  algorithm  implementation.  One  could  come  up  with  some  clever  scheme  to  predict  or  esti¬ 
mate  the  color  ratio,  but  it  will  be  shown  that  the  performance  of  the  Max_SNR  is  not  sufficiently 
better  than  the  other  algorithms  to  warrant  the  effort.  The  other  algorithms  do  not  require  this  infor¬ 
mation,  and  are  much  simpler  to  implement.  The  Max_SNR  algorithm  serves  as  an  upper  bound  to 
the  performance  obtainable  from  a  weighted-difference  algorithm. 

The  weight  for  variant  2,  <umin_noise,  can  be  determined  by  taking  the  derivative  of  that  quantity 
and  finding  the  value  of  co  that  minimizes  the  denominator.  There  is  only  one  solution  for  co  in  this 
case  which  is 


C^min_noise  =  P  •  (12) 

The  Min_Diff  algorithm  does  not  require  knowledge  of  the  target  color.  However,  the  algorithm  has  a 
similar  problem  to  Max_SNR  in  that  it  must  know  which  band  has  the  largest  input  SNR.  If  this  is 
unknown,  the  algorithm  must  be  performed  twice,  once  for  x  hypothesized  as  the  largest  and  once  for 
y.  Both  outputs  must  be  threshold  tested  for  detection.,.  The  penalties  for  not  knowing  which  channel 
has  the  largest  input  SNR  are  increased  computation  and  either  a  higher  false-alarm  rate  or  a  higher 
threshold  value— which  will  reduce  the  number  of  detections. 

The  weight  for  the  simple  difference  variant  is  simply 

tWsim_di((  =  1*0  •  (13) 

The  output  SNR0Ut  for  each  variant  can  be  found  by  substituting  the  appropriate  solutions  for  co  back 
into  equation  (8).  Table  1  summarizes  the  three  variants  of  the  weighted-difference  algorithm. 
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Table  1.  Variants  of  the  weighted-difference  algorithm. 


Casel :  Max_SNR 


SNR 


S?  .  (1  - 

o?  1  -  2(o  P+  a)2 


■» 


Case  2:  Min  Noise 


(0  =  p 


p2 

SNR  =  -3 
Ol 


(1  -PC)2 
1  -P2 


Case  3:  Sim  Diff 


(0  =  1.0 


c2 

SNR  =  -j 

Oi 


(1  -  c)2 

2(1  -p) 


Figure  2  is  a  plot  of  the  theoretical  output  SNR  for  a  band-to-band  correlation  coefficient  of 
p  =  0.9.  The  dB  gain  is  very  dependent  on  the  color  ratio.  When  the  color  ratio  is  1,  there  is  no 
color  to  the  target.  When  the  color  ratio  is  negative,  the  target  contrasts  in  the  two  bands  have  oppo¬ 
site  signs;  i.e.,  the  target  in  one  band  is  above  its  mean  level  and  the  target  in  the  other  band  is 
below  its  mean  level.  The  gain  from  these  multispectral  algorithms  is  greatest  for  negative  color  ratios. 
This  is  an  important  point  since,  for  this  to  occur,  the  target  must  be  near  zero  contrast  and  difficult 
to  detect  by  spatial  processing  alone. 


SNR  dB  GAIN  FOR  SPATIAL  DIFF 
BAND-TO-BAND  CORREL  COEFF  0.9 


COLOR  RATIO 

-  MAX_SNR  -  MIN_NOISE  -  SIM_DIFF 

Figure  2.  Weighted-difference  algorithm  variants  for  a  background 
band-to-band  correlation  of  0.9. 
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Figure  3  shows  three  cases  of  interest.  The  first  shows  the  radiance  levels  when  the  target  appears 
hot  compared  to  the  background  mean  level.  In  this  case,  the  input  SNRs  are  positive  and  large  (for 
this  simple  argument,  the  variance  of  the  background  about  the  mean  level  is  not  shown) .  The  SNRs 
in  the  two  channels  are  different  only  by  the  factor  of  the  emissivity.  Thus  the  ratio,  c,  will  be  positive 
and  close  to  1.  In  case  2,  the  input  signal-to-noise  ratios  are  small  because  they  are  close  to  the  back¬ 
ground  mean  level.  Here  the  ratio  of  the  SNRs  can  be  magnified,  depending  on  the  mean  level.  The 
ratio  of  the  SNRs  could  be  large  and  possibly  negative.  This  condition  results  in  large  SNR  gain.  In 
case  3,  the  target  appears  cold  compared  to  the  background.  If  the  SNRs  are  large  (but  negative),  the 
ratio  again  will  be  positive  and  close  to  1.  This  results  in  a  small  gain. 

Figure  4  is  a  theoretically  derived  curve  comparing  the  output  SNR  from  a  spatial-differenced  sys¬ 
tem  with  the  best  input  SNR.  The  model  parameters  were  arbitrarily  chosen,  but  the  figure  illustrates 
the  behavior  of  the  differenced  algorithm.  It  tends  to  provide  gain  when  the  input  SNR  is  small.  This 
characteristic  “iills  in”  the  SNR  output  curve  when  the  input  SNR  is  small. 


3.0  COMPARISON  OF  THEORETICAL  AND 
EXPERIMENTAL  RESULTS 

In  this  section,  the  theoretical  performance  is  compared  to  measured  performance  for  two  multi- 
spectral  infrared  scenes  collected  by  the  NASA  TIMS  sensor.  The  first  scene  was  collected  over  a 
rural  area  near  Adelaide,  Australia,  and  the  second  scene  was  collected  over  a  mixed  wooded  and 
open-field  area  near  the  Redstone  Arsenal  in  Huntsville,  Alabama.  The  TIMS  sensor  has  six  channels 
between  8  and  12  urn;  however,  only  two  were  used  for  weighted-difference  application.  Scenes  were 
selected  that  were  512  x  512  and  256  X  256  pixels  in  dimension  for  the  Adelaide  scene  and 
256  X  256  for  the  Redstone  scene.  Each  scene  contains  a  number  of  manmade  objects  in  a  cluttered 
background  that  served  as  targets.  There  were  52  objects  (houses,  farm  buildings,  watertanks,  etc.)  in 
the  Adelaide  512  x  512  scene  that  were  declared  targets  and  used  to  evaluate  the  dB  gains  as  a 
function  of  the  input  SNR  and  detection  statistics.  There  were  14  declared  targets  in  the  Adelaide  256 
X  256  image  that  were  used  to  evaluate  SNR  statistics.  The  analysis  was  repeated  on  the  Redstone 
scene;  however,  this  time  military  tanks  (M60  and  M48)  were  located  in  the  scene  as  part  of  the  ex¬ 
periment.  The  Redstone  scene  was  used  to  confirm  the  results  achieved  by  using  the  Adelaide  data. 


3.1  EXPERIMENTAL  RESULTS  FROM  THE  ADELAIDE  SCENE 

The  TIMS  Australian  data  set  was  acquired  as  part  of  a  joint  U.S  and  Australian  measurements 
program  conducted  in  the  fall  of  1985  over  several  regions  of  Australia.  The  data  set  was  obtained 
from  a  C-130  aircraft  at  a  variety  of  altitudes.  A  review  of  the  data  taken  during  this  particular  acqui¬ 
sition  found  a  rural  natural  background  area  in  the  hills  east  of  Adelaide.  A  subsection  of  this  area 
containing  structures  was  selected  for  spectral  processing  [6]. 

Figure  5  is  a  grey-level  image  of  the  Adelaide  256  x  256  scene.  The  C-130  flew  at  a  low  altitude, 
which  resulted  in  a  ground  resolution  of  2.5-meter  pixels.  Channels  1  and  5  were  selected  for  process¬ 
ing  because  they  gave  the  best  target  color  ratios  and  band-to-band  clutter  correlation.  Figure  6  shows 
the  output  scene  from  the  weighted-difference  algorithm.  The  minimum  noise  variant  with  full 
preprocessing  (i.e.,  line-by-line  normalization  and  high-pass  filtering)  was  used.  Note  that  the  output 
targets  can  have  positive  (bright)  or  negative  (dark)  contrast.  Details  of  the  numeric  results  are  given 
in  appendix  A. 
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SNR  (dB)  =  10  logtpsnr(7)]2 
Let  es\-  0.80.  effective  emissivity  of  target  in  band  1 
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Band'to-band  correlation  coefficient  p=  0.9. 

Figure  4.  Channel  radiance  versus  temperature. 
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A  rnrai  area  near  Adelaide,  South  Australia,  consisting  of  scattered  trees,  a  forested 
region,  rolling  hills,  and  manmade  structures.  The  image  has  256  X  256  pixels  that 
have  2.5-meter  ground  resolution  Channels  1  and  5  were  used  for  multispectral  proc¬ 
essing. 


Figure  5.  Adelaide  scene. 
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The  results  of  processing  the  Adelaide  256  x  256  data  are  compared  to  the  theoretical  output 
signal-to-noise  ratio  in  figures  7a  and  7b.  First  the  muhispectral  images  were  processed  without 
preprocessing  and  the  results  are  shown  in  figure  7a.  The  signal-to-noise  gain,  measured  in  dB,  is  plot¬ 
ted  as  a  function  of  the  color  ratio.  The  correlation  coefficient  for  the  Adelaide  scene  was  0.9929. 

* 

There  are  three  curves  on  the  graph  that  are  the  theoretically  predicted  gain  for  the  three  variants  of 
the  weighted-difference  algorithm.  The  symbols  represent  measured  gain  on  targets  in  the  Adelaide 
scene-  There  is  excellent  agreement  between  the  theory  and  measured  results. 

Similarly,  figure  7B  presents  the  same  information,  except  that  in  this  case  the  data  set  was 
preprocessed  as  shown  in  figure  1.  The  signal-to-noise  ratio  gain  is  measured  from  the  input  to  the 
output  of  the  weighted  difference  algorithm,  and  does  not  include  the  SNR  gain  or  loss  caused  by  the 
preprocessing.  In  general,  the  weighted  difference  gains  are  much  lower  with  preprocessing  than  with¬ 
out  preprocessing,  as  can  be  seen  by  comparing  the  results  to  the  example  in  figure  7.  The  reason  for 
the  drop  in  gain  is  that  the  preprocessing  has  reduced  both  the  band-to-band  correlation  and  the 
color  ratio.  However,  the  preprocessing  also  provides  gain,  so  the  overall  effect  is  a  net  improvement 
over  no  preprocessing.  The  output  signal-to-noise  ratio  level  for  preprocessing  and  no  preprocessing 
are  summarized  in  tables  A4  through  A6  of  appendix  A.  With  a  few  exceptions,  the  combination  of 
preprocessing  and  weighted  difference  improved  the  overall  output  signal-to-noise  ratio  by  1  or  2  dB. 
The  significant  conclusion  from  figures  7a  and  7b  is  that  the  SNR  performance  of  the  weighted- 
difference  algorithm  can  be  predicted. 

Figure  8  shows  the  output  signal-to-noise  ratio  plotted  versus  the  input  average  pixel-to-noise  ratio. 
The  pixel  signal-to-noise  ratio  is  given  by  equation  (4).  The  average  psnr  is  a  function  of  the  differ¬ 
ence  between  the  target  and  the  mean  background  temperatures.  Using  measured  data  from  the 
Adelaide  512  X  512  scene  (see  appendix  C),  figure  8  compares  the  best  input  signal-to-noise  ratio  to 
the  min_noise  (without  preprocessing)  output  SNR.  The  data  have  the  same  trend  as  the  theoretical 
data  shown  in  figure  4.  The  data  are  not  smooth  because  the  targets  do  not  have  the  same  emissivity, 
which  causes  the  gain  to  vary  for  targets  even  at  the  same  temperature.  Regardless,  the  gain  appears 
to  be  maximum  near  zero,  as  in  figure  4,  and  shows  the  characteristic  of  the  weighted-difference  algo¬ 
rithm  to  “fill  in”  the  output  SNR  when  the  target  is  small  because  of  marginal  target  temperature  dif¬ 
ference  from  the  background. 

The  detection  statistics  for  the  52  targets  in  the  Adelaide  512  x  512  scene  were  computed  for 
10"3  and  10~4  false-alarm  rates,  and  the  results  are  summarized  in  table  2.  Details  of  the  numeric 
results  are  provided  in  appendix  C.  The  threshold  for  the  false-alarm  rates  was  computed  from  the 
data  after  the  52  targets  were  omitted  from  the  512  x  512  image.  As  shown  in  figure  8,  the  perform¬ 
ance  improvement  of  the  weighted-difference  algorithm  depends  on  the  relative  target  temperature  to 
the  average  background  temperature  as  well  as  the  target  emissivity.  Thus  the  number  of  detections  in 
table  2  depend  on  the  statistical  distribution  of  the  target  temperature  and  emissivity  provided  in  the 
Adelaide  scene. 


Table  2.  Detection  statistics  for  52  targets  in  the  Adelaide  scene. 


Processing  Algorithm  Detections  at  CFAR 

(Bands  1  and  5)  10-3  1(M 


No  Preprocessing 

Band  5— Raw  Input .  28  20 

Simple  Difference  .  44  18 

Minimum  Noise .  45  18 

Preprocessed 

Band  5 .  34  19 

Simple  Difference  .  47  23 

Minimum  Noise .  ^7  19. 
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BAM)-10-BAND  CORRELCOEF  0.9929 


0.30  0.50  0.70  0.90 

COLOR  RATIO 


o  MAXjSNR  X  MIN_NOISE  SIM_DIFF 


(a) 


BAND-TO-BAND  CORREL  COEF  0.988112 


COLOR  RATIO 


o  MAXSNR  X  MIN_NOISE  *  SIM_DIFF 


<b) 


Actual  data  points  from  the  Adelaide  256  X  256  scene  plotted  on  the  theoretical  curves 
for  the  three  variants  of  the  weighted-difference  algorithm. 


Figure  7a-b.  dB  gain  as  a  function  of  color  ratio. 
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Actual  data  points  from  the  Adelaide  512  X  512  scene.  Compare  with  figure  4. 


Figure  8.  dB  gain  as  a  function  of  input  SNR. 
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Without  any  processing,  28  of  the  52  targets  could  be  detected  at  a  1CT3  false-alarm  rate.  With 
lull  preprocessing  and  the  minimum-noise-difference  algorithm,  this  number  was  increased  to  47.  The 
detection  rate  improvement  from  54%  to  90%  is  a  very  significant  result. 

3.2  EXPERIMENTAL  RESULTS  FROM  THE  REDSTONE  SCENE 

The  Redstone  experiment  was  performed  on  9  March,  1989,  at  Redstone  Arsenal,  Alabama,  using 
the  TIMS  sensor.  A  Lear  Jet  was  flown  during  two  time  periods,  08C0  to  0900  and  1100  to  1200  h. 

A  data  set  was  taken  during  each  of  these  flights  at  three  different  altitudes,  resulting  in  a  resolution 
of  2.5,  5,  and  10  meters.  The  weather  was  clear,  with  a  slight  breeze  and  moderate  haze  that  burned 
off  during  the  morning.  The  air  temperature  increased  from  39  to  51°F, 

The  experimental  setup  was  located  at  Redstone  Arsenal’s  Test  Area  No.  3,  which  consists  of 
slightly  rolling  hills  with  a  grass-covered  stretch  of  land  surrounded  by  loblolly  pines  and  deciduous 
trees.  The  grass  was  disturbed  as  a  result  of  tank  movements  in  the  field.  Three  tanks  were  set  up  as 
targets,  and  nearby  temporary  buildings  and  a  concrete  pad  located  to  the  south  of  the  tanks  were 
used  as  targets  of  opportunity.  Tank  A,  an  M48,  was  located  on  the  edge  of  the  field  with  trees  just 
to  the  west.  Tank  B,  an  M60,  was  under  the  trees  on  a  dirt  road  covered  with  pine  needles.  The 
trees  provided  about  a  50%  canopy  over  the  tank.  Tank  C,  another  M48,  was  located  in  the  open 
field.  The  ground  was  less  disturbed  around  tank  C  than  tank  A.  Tank  B,  the  tank  under  the  trees, 
was  not  observed  with  the  infrared  sensor  or  high-resolution  photography. 

Both  tank  A  and  tank  C  were  detected  in  the  2.5-m  data  by  means  of  several  algorithms  in  both 
the  morning  and  noon  scenes.  Tank  C,  the  open-area  tank,  was  also  detected  in  the  lower-resolution 
data.  In  this  section,  we  will  concentrate  on  the  detections  made  during  the  early-morning  2.5-m 
scene,  with  particular  emphasis  on  the  tree-line  tank  A.  A  gray-level  image  of  the  scene  showing  the 
tank  locations  is  given  in  figure  9. 

Bands  1  and  4  of  the  TIMS  data  were  spectrally  processed  and  the  detailed  numeric  results  are 
given  in  appendix  D.  Bands  1  and  4  were  chosen  since  they  show  spectral  color.  The  results  for  a 
32  x  32  image  centered  on  the  tree  line  of  tank  A  are  shown  in  figure  10.  The  tree-line  tank,  not 
discernible  in  the  input  data  in  figure  10A,  can  be  seen  easily  in  figure  10B. 

The  location  of  a  target  detection  relative  to  false  alarms  is  an  important  indication  of  algorithm 
performance,  since  a  constant  false-alarm  rate  system  would  need  to  suppress  false  alarms  to  a  speci¬ 
fied  level.  The  cumulative  distributions  for  the  unprocessed  input  and  processed  data  are  shown  in 
figures  11a  and  lib.  Tank  A  and  C  pixels  were  removed  before  calculating  the  distributions.  The 
curves  can  be  used  to  estimate  the  false-alarm  rate  that  would  result  for  a  given  threshold  decision 
level.  For  unprocessed  input  data,  the  false-alarm  rate  would  have  to  be  2  X  10"1  before  the  target 
could  be  detected.  After  processing,  a  more  reasonable  false-alarm  rate  of  only  7  X  10'5  would  have 
to  be  tolerated.  (Recall  that,  if  the  strongest  input  channel  were  not  known  a  priori,  a  double-sided 
threshold  would  be  required).  The  double-sided  threshold  could  double  the  number  of  false  alarms  to 
1.4  x  10‘4,  which  is  still  a  reasonable  value.  The  false  alarms  at  thresholds  greater  than  tank  A  were 
caused  by  man-made  structures  (concrete  pad,  buildings,  vehicles)  in  the  scene,  which  were  also  col¬ 
ored.  These  man-made  objects  can  be  considered  as  targets  in  the  context  of  a  natural  scene. 

The  various  colored,  manmade  objects  in  the  scene  can  be  used  to  verify  the  theory  presented 
earlier.  The  output  depends  on  the  amount  of  target  coloring  [see  equations  (8)  and  (10)].  Each  of 
these  objects  has  a  different  color  ratio.  These  objects  were  identified  by  means  of  site  visits  and  aer¬ 
ial  photography.  The  objects  investigated  included  two  structures,  two  pixels  on  tank  C,  tank  A,  and  a 
rectangular  concrete  pad. 
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Sceno  consists  of  slightly  rolling  hills,  with  a  grass-covered  stretch  of  land 
surrounded  by  loblolly  pines  and  deciduous  trees. 


Figure  9.  Huntsville,  Alabama— Redstone  scene. 
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These  results  are  presented  for  the  three  algorithm  variants  in  figure  12.  The  solid  lines  of  figure 
12  are  plots  of  the  theoretical  output  SNR  for  a  band-to*band  correlation  coefficient  of  p  =  0:9.  The 
dB  gain  is  very  dependent  on  the  color  ratio.  When  the  color  ratio  is  1,  there  is  no  color  to  the  tar¬ 
get.  When  the  color  ratio  is  negative,  the  target  contrasts  in  the  two  bands  have  opposite  signs;  i.e., 
the  target  in  one  band  is  above  its  mean  level  and  the  target  in  the  other  band  is  below  its  mean 
level. 

The  symbols  represent  data  points  taken  from  targets  in  the  Redstone  scene.  For  these  data 
points,  the  gain  is  plotted  against  the  input  color  ratio.  The  data  set  shows  an  excellent  match  to  the 
theory. 


SNR  dB  GAIN  FOR  SPATIAL  DIFF 
BAND-TO-BAND  CORREL  COEFF  O.S 


Actual  data  points  from  the  Redstone  scene  plotted  on  the  theoretical  curves  for  the  three  variants 
of  the  weighted-difference  algorithm.  For  this  data  set,  preprocessing  was  not  performed. 

Figure  12.  dB  gain  as  a  function  of  color  ratio. 


4,0  CONCLUSIONS 

In  this  report,  the  weighted-difference  algorithm  was  derived,  analyzed,  and  compared  to  meas¬ 
ured  data.  Three  variants  of  the  algorithm  were  evaluated:  (1)  simple  difference;  (2)  minimum  noise; 
and  (3)  maximum  output  SNR.  The  variants  require  increasing  amounts  of  prior  information  to  be 
implemented.  Although  the  corresponding  performance  of  the  variants  also  increases,  the  increase  is 
very  small  for  targets  that  are  in  the  detectable  range  of  input  SNR  and  color  ratio. 

The  performance  of  the  weighted-difference  algorithm  depends  on  the  correlation  coefficient 
between  the  two  spectral  bands,  the  ratio  of  the  target  pixels  (color  ratio),  and  the  input  target’s 
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signal-to-noise  ratio.  In  this  report,  the  theoretical  performance  has  been  compared  to  measured  per¬ 
formance  for  two  scenes.  The  scenes  were  collected  by  the  NASA  TIMS  sensor  over  a  rural  area  near 
Adelaide,  Australia  and  over  a  wooded  area  near  the  Redstone  Arsenal.  Scenes  were  selected  that 
were  512  X  512  and  256  X  256  pixels  in  dimension  and  contained  a  number  of  targets  in  a  cluttered 
background.  There  were  52  objects  (buildings,  watertanks,  etc.)  in  the  Adelaide  512  x  512  scene 
that  were  declared  targets  and  used  to  evaluate  the  dB  gain  as  a  function  of  input  SNR  and  detection 
statistics.  There  were  14  declared  targets  in  the  Adelaide  256  X  256  image  that  were  used  to  evaluate 
SNR  statistics.  The  analysis  was  repeated  on  the  Redstone  scene;  however,  this  time  there  were  mili¬ 
tary  tanks  (M60  and  M48)  in  the  scene.  The  Redstone  scene  was  used  to  confirm  the  results 
achieved  by  using  the  Adelaide  scene. 

The  theoretical  and  measured  results  agree  extremely  well.  For  a  given  correlation  coefficient  and 
color  ratio,  the  amount  of  signal-to-noise  ratio  gain  can  be  predicted.  However,  target  input  SNRs  and 
color  ratios  can  vary  considerably.  For  the  targets  and  scenes  evaluated  here,  the  typical  gains 
achieved  ranged  from  a  few  dB  loss  (targets  without  color)  to  a  maximum  of  approximately  20  dB. 
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APPENDIX  A:  DETAILS  OF  NUMERIC  RESULTS, 
ADELAIDE  DATA  SET 


The  TIMS  Adelaide  data  set  was  processed  to  compare  the  results  to  the  spatial  difference  model 
developed  by  L.  E.  Hoff.  The  first  step  in  this  processing  was  to  determine  which  two  of  the  six 
TIMS  bands  met  the  criterion  of  high  correlation  between  the  bands.  A  routine  was  run  on  the  raw 
data  that  produces  a  correlation  coefficient  matrix  as  shown  in  figure  A-l.  It  was  also  desirable  for 
this  processing  that  the  bands  be  widely  separated,  which  tends  to  preserve  target  coloring.  Bands  1 
and  5,  with  a  correlation  coefficient  of  0.9929,  were  chosen. 
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Figure  A-l.  Correlation  coefficient  matrix. 


The  second  step  in  the  processing  was  to  window  the  area  of  interest.  An  area  that  was 
256  X  256  in  size  was  chosen.  From  this  area,  14  target  pixels  were  chosen  for  this  processing  run. 
Figure  A-2  identifies  the  target  pixels  by  line  and  column. 
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Figure  A-2.  Adelaide  target  pixels. 


Figure  A-3  outlines  the  processing  steps  used  on  the  256  X  256  image.  This  diagram  shows  the 
general  data  flow  but  does  not  show  the  switches  that  can  be  made  in  preprocessing  of  the  data. 
Three  different  preprocessing  methods  were  used  on  this  data  set.  The  first  case  did  a  global 
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normalization  on  the  image,  which  took  the  mean  and  variance  of  the  entire  256  X  256  image  to 
compute  pixel  values.  The  second  case  did  a  line-by-line  normalizadon  on  the  data,  which  entailed 
taking  the  mean  and  variance  of  each  line  to  determine  the  pixel  values  of  each,  and  then  doing  a 
global  normalizadon  on  the  data.  The  third  case  did  a  line-by-line  normalizadon  on  the  data,  followed 
by  a  high-pass,  5x5  filter  on  the  data.  As  a  last  step,  the  data  set  was  globally  normalized.  The 
final  output  of  the  preprocessing  phase  of  this  image  processing  consisted  of  six  files,  three  each  of 
band  1  and  band  5,  globally  normalized,  line-by-line  normalized,  and  filtered. 


I 


■A 


RAW  PREPROCESSING  WEIGHTED-DIFFERENCE  ALGORITHM  PROCESSED 

DATA  DATA 


Figure  A-3.  Processing  steps  used  on  256  x  256  image. 


For  each  case,  and  at  each  target  pixel,  the  band  with  the  largest  signal-to-noise  (SNR)  was 
determined  and  this  data  set  was  entered  into  a  database  (see  tables  A-l,  A-2,  and  A-3).  The  pixel 
values  that  determined  the  color  ratio  for  comparison  purposes  were  also  entered  at  this  time. 

On  completion  of  the  preprocessing  steps,  there  were  three  weighted-difference  algorithms  applied 
to  each  set  (band  1  and  band  5)  of  data.  Each  of  these  algorithms  used  a  different  weighting  factor; 
one  maximized  the  output  SNR,  one  minimized  the  noise,  and  one  used  a  weighting  factor  of  1. 
These  algorithms  are  referred  to  as  Max_SNR,  Min_Noise,  and  Sim_Diff  (simple  difference).  The 
Min_Noise  algorithm  was  applied  twice  to  each  data  set.  One  run  applied  band  1  to  band  5,  and  the 
second  run  applied  band  5  to  band  1.  The  SNR  for  each  target  pixel  from  each  output  image  was 
then  measured  and  entered  into  the  database  with  the  input  band’s  SNRs  (tables  A-l  through  A-3). 
The  dB  gain  was  then  computed  from  the  input  and  output  SNR  information  and  entered  into  the 
database.  The  data  provided  the  necessary  information  for  comparing  the  spatial-difference  model  to 
the  Adelaide  data. 

Tables  A-4,  A-5,  and  A-6  consolidate  the  SNR  data  for  each  weighted  difference  algorithm  for 
quick-look  purposes. 
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Table  A-l.  Adelaide  normalizsd  (entire  image  only)  target  data. 
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Table  A-2.  Adelaide  normalized  (line  by  line,  then  entire  image)  target  data. 
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Table  A-3.  Adelaide  (normalized  line  by  line/high-pass  filter/then  renormalized)  target  data. 
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Table  A-4.  SNRs  for  various  processing  (normalized  only,  normalized  line 
by  line,  filtered)  maximum  SNR— Adelaide  channels  1  and  5  (8/15/90). 


>  SNR  IN  (dB) 

NORMALIZED 

NORMALIZED 

HIGH  PASS 

NAME 

LINE 

PIXEL 

BAND 

VALUE 

WHOLE  IMAGE 

LN  X  LN 

FILTER 

ala 

133 

50 

5 

13.89994 

19.00481 

17.48344 

20.51602 

alb 

133 

51 

5 

14.12785 

22.89334 

21-98891 

24.43176 

a2a 

137 

76 

5 

14.12785 

22.22293 

19.72355 

20.87023 

a2b 

137 

77 

5 

17.7531 

23.41788 

20.28915 

22.44913 

a3a 

162 

51 

5 

2.05 

13.45629 

14.23086 

15.91557 

a3b 

163 

50 

5 

4.4093 

13.59113 

13.36969 

15.76314 

a4a 

57 

194 

5 

15.76334 

21.0701 

20.88108 

21.90009 

a4b 

58 

195 

5 

16.47845 

23.18069 

23.11026 

24.32445 

a5a 

237 

66 

5 

14.12785 

22.22293 

21-68434 

23.67152 

a5b 

237 

67 

5 

10.895 

18.45592 

17.93947 

19.87536 

a6a 

245 

158 

5 

4.92731 

20.73613 

19.63205 

18.62448 

a6b 

246 

158  | 

5 

10.13666 

21.66016 

20.81269 

20.45167 

a7 

199 

92 

1 

14.41609 

15.29576 

14.22003 

15.38112 

a8 

188 

235 

5 

10.82128 

20.88979 

20.65253 

21.27448 

Table  A-5.  SNRs  for  various  processing  (normalized  only,  normalized  line 
by  line,  filtered)  minimum  noise— Adelaide  channels  1  and  5  (8/15/90). 


NAME 

LINE 

PIXEL! 

>  SI 
BAND 

$R  IN  (dB) 
VALUE 

NORMALIZED 
WHOLE  IMAGE 

NORMALIZED  I 
LN  X  LN 

HIGH  PASS 
FILTER 

ala 

133 

50 

5 

13.89994 

17.39955 

15.73921 

17.98848 

alb  [ 

133 

51 

5 

14.12785 

22.27282 

21.41471 

23.46915 

a2a 

137 

76 

5 

14.12785 

21.49039 

18.90941 

18.97939 

a2b 

137 

77 

5 

17.7531 

22.04127 

18.38488 

18.04218 

a3a 

162 

51 

5 

2.05 

13.12955 

14.04098 

14.86855 

a3b 

163 

50 

5 

4.4093 

13.03079 

12.94531 

13.83025 

a4a 

57 

194 

5 

15.76334 

19.55179 

19.23371 

18.99117 

a4b 

58 

195 

5 

16.47845 

22,13463 

21.9855 

22.44597 

a5a 

237 

66 

5 

14.12785 

21,49039 

21.01014 

22.34047 

a5b 

237 

67 

5 

10.895 

17.61704 

17.15941 

18.30028 

a6a 

245 

158 

5 

4.92731 

20.62004 

19.56133 

18.57658 

a6b 

246 

158 

5 

10.13666 

21.3422 

20.58337 

20.12804 

a7 

199 

92 

1 

14.41609 

9.82421 

8.52137 

12.93983 

a8 

188 

235 

5 

10.82128 

20.43851 

20.28568 

20.96671 
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Table  A-6.  SNRs  for  various  processing  (normalized  only,  normalized  line  by 
line,  filtered)  simple  difference— Adelaide  channels  1  and  5  (8/15/90). 


I >  SNR  IK  (dB)  NORMALIZED  I  NORMALIZED 

NAME  LINE  PIXEL  BAND  VALUE  WHOLE  IMAGE  LN  X  LN 

ala  133  50  5  13.89994  17.72324  16.05657 

alb  133  51  5  14.12785  22.45768  21.57983 

a2a  137  76  5  14.12785  21.69389  19.11112 

a2b  137  77  5  17.7531  22.33606  18.72013 

a3a  162  51  5  2.05  13.25748  14.12847 

a3b  163  50  5  4.4093  13.20508  13.08428 

a4a  57  194  5  15.76334  19.8646  19.54006 

a4b  58  195  5  ;  16.47845  22.38491  22.22903 

a5a  237  66  5  14.12785  21.69389  21.19112 

a5b  237  67  5  10.895  17.83706  17.35621 

a6a  245  158  5  4.92731  20.68923  19.60928 

a6b  246  158  5  10.13666  21.46809  20.68101 

a7  199  92  1  14.41609  8.94525  6.68011 

a8  188  235  5  10.82128  20.59246  20.41357 


HIGH  PASS 
FILTER 
18.5398 
23.77154 
19.43594 
18.86422 
15.18646 
14.29322 
19.59766 
22.90062 
22.70777 
18.70749 
18.62112 
20.28698 
12.31004 
21.12094 
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APPENDIX  B:  MODEL  FOR  TWO-CHANNEL  SYSTEM 


The  following  is  a  model  for  a  two-channel  system.  Let 

Rk  =  Sk  or  Nk  k  =  l,2 

be  the  received  radiance  in  band  k.  If  signal  is  present,  the  background,  TV* ,  is  occluded.  The  spatial 
difference  algorithm  is  formulated  as 


For  the  signal  pixel  (s),  we  get 


d(R )  —  Ri  —  0)R2  . 


d(S)  =  Si  -  wS2 

and  for  the  background  pixels  (noise  and  clutter),  we  get 

d(N)  =  N%  +  coN2  . 
The  output  SNR  is  defined  as  follows 


\m  1 

var[d(A0] 


Assume  that  each  channel  input  has  been  normalized  so  that  image  mean  is  zero  and  the  variance  is 
1.0.  Define  the  following  two  parameters: 

S2 

c  =  —  ,  where  S2<S\  ,  -  1  <  c  <  1 

S 1 


Cov  [N\,N2] 

Oi02 


First,  let  us  derive  the  SNR: 

var[d (W) ]  =  E[(N\  -  a>N2)2] 

=  E[N\  -  2a)N\N2  +  o?N\] 
=  of  -  2(0fXJiO2  +  to2of 

=  1  -  2ci)  p  +  a>2 
and 


giving 


Sni  - 


(Si  ~  m,) 

Oi 


i  =  1,2 


SNR 


(Sni  ~  0)Sni)2 
(1  -  2a)p  +  ft)2) 


•&(!  -  0)C)2  _ 

1  -  2ft)  p  +  ft)2 
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We  derive  a>  three  ways:  (1)  choose  cd  to  maximize  SNR;  (2)  choose  o>  to  minimize  V  ar[d(N)); 
(3)  choose  o)  equal  to  1.  To  find  the  cd  that  maximizes  the  SNR,  we  perform  the  following  deriva¬ 
tions: 


aSNR 

do) 


2&(1  -  oc)(-c)  #i(l  -  cdc)2 
1  -  2(d  p  +  a)2  (1  -  2a)  p  +  ci>2)2 


2p)(-l) 


-■Snl(l  -  0>c)  j  2c  +  (1  -  CPC) (2ft)  -  2p) 
1  -  2cd  p  +  cd2  |_  1-  2a>p  +  cd2 

~ 2Sni(l  -  (oc)  |”c  +  (1  -  cdc) (ft)  -  p)~l 
1  -  2a)p  +  cd2  |_  1  -  2cip  +  CD2  J  * 


The  above  equation  has  two  zeros.  The  first  one,  easily  seen  in  the  first  term,  is  cd  =  1/c.  The  sec¬ 
ond  one  is  derived  as  follows: 

(1  -  cd c)(o>  -  p)  a)  -  a)2c  -  p  +  ocp  +  c  -  2q>cp  +  cd2c 

1  -  2a)  p  +  cd2  1  -  2cdp  +  CD2 


CD  +  C  -  (OC  P  -  p 

1  -  2cdp  +  a >2 


0)(1  -  cp)  +  C  -  p 
1  -  2cdP  +  CD2 


which,  when  we  set  equal  to  zero,  gives 


0) 


P-  c 
1  -  pc 


The  first  term  gives  a  minimum,  the  second  a  maximum.  Figure  B-l  is  a  plot  of  SNR  as  a  function  of 
o),  showing  the  maximum  and  minimum.  This  plot  has  p  equal  to  0.99  and  c  equal  to  0.9.  For  case  2 
we  take  the  derivative  of  the  noise  term  with  respect  to  (o . 

dVar[d(N )] 

- ‘  =  -  2  p  -  2  CD  =»p=CD 

do) 


and 


SNR  -  Snx 


(1  -  PC)2 
1  -p  2 


For  case  3,  we  just  set  cd  =  1,  giving 


SNR 


&c(l  -  c)2 

2(1  -  p)  • 


i 
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APPENDIX  C:  ANALYSIS  OF  DATA,  ADELAIDE 
512  x  512  SCENE 


The  TIMS  Adelaide  512  x  512  data  set  was  analyzed  to  evaluate  detection  statistics  for  10'3  and 
10‘4  false-alarm  rates  as  a  follow-on  to  investigation  of  weighted-difference  signal-processing  algorithms 
for  detecting  g  und  targets  by  using  dual-band  IR  data.  Extensive  investigation  of  these  algorithms 
was  done  by  using  the  TIMS  Adelaide  256  X  256  data  set;  however,  for  determining  detection  statis¬ 
tics,  the  TIMS  Adelaide  512  X  512  image  was  chosen  because  it  contained  significantly  more  targets. 
Figure  C-l  identifies  these  targets  by  row  and  column  number. 
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Figure  C-l.  TIMS  Adelaide  512  x  512  image  targets. 


The  detection  statistics  were  determined  for  two  sets  of  data:  data  that  had  not  been  preprocessed 
and  preprocessed  data.  The  preprocessed  data  set  was  line-by-line  normalized,  filtered  with  a  high- 
pass  filter,  and  then  globally  normalized.  There  were  two  weighted-difference  algorithms,  simple  differ¬ 
ence  and  minimum  noise,  applied  to  each  set. 

Histograms  were  formed  for  each  image  after  omitting  the  52  targets.  Two  accumulated  statistics 
were  obtained  for  each  histogram:  ascending  values  and  descending  values.  These  accumulated  statis¬ 
tics  were  used  to  estimate  the  false-alarm  rate  associated  with  each  target.  The  threshold  for  the  false- 
alarm  rates  was  then  obtained  from  the  data.  Attachment  (1),  which  follows  table  C-10,  contains  the 
numeric  results  of  this  computation.  At  the  predetermined  thresholds  of  10"3  and  10*4  ,  the  number 
of  detections  that  occurred  was  determined  by  using  previously  measured  pixel  values  in  tables  C-l 
through  C-10.  This  determination  included  targets  that  were  positive,  with  a  pixel  value  larger  than  the 
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image  mean,  and  negative  targets,  with  a  pixel  value  less  than  the  mean.  Figure  C-2  summarizes  the 
results. 


Processing  Algorithm 
(Bands  1  and  5) 

Detections  at  CFAR 
1<H  10-* 

No  Preprocessing 

Band  5— Raw  Input . 

28 

20 

Simple  Difference  . 

44 

18 

Minimum  Noise . 

45 

18 

Preprocessed 

Band  5 . 

34 

d9 

Simple  Difference . 

47 

23 

Minimum  Noise . 

47 

19 

Figure  C-2.  Detection  statistics  for  52  targets  in  the  Adelaide  scene. 
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Table  C-l.  Channel  1,  no  preprocessing. 


Line 

Name 

Pixel  val 

—  ill  1  III  1  1 

133 

50 

A1A 

59.00000 

11.64553 

10.83174 

133 

51 

A1B 

63.00000 

10.68036 

9.76055 

137 

76 

A2A 

62.00000 

12.22594 

10.04098 

137 

77 

A2B 

43.00000 

16.63482 

14.14314 

162 

51 

A3A 

84.00000 

-1.43988 

-0.38641 

163 

50 

A3B 

81.00000 

1.53535 

2.00270 

57 

194 

A4A 

52.00000 

12.56733 

12.43736 

58 

195 

A4B 

51.00000 

12.63880 

12.64430 

237 

66 

A5A 

62.00000 

10.32106 

10.04098 

237 

67 

A5B 

70.00000 

7.74971 

7.49383 

245 

158 

A6A  . 

93.00000 

-23.70481 

-26.37514 

246 

158 

A6B 

103.00000 

3.53570 

-0.34164 

199 

92 

A7 

129.00000 

10.89036 

11.09208 

188 

235 

A8 

106.00000 

2.80526 

2.03672 

259 

119 

A9 

64.00000 

11.30198 

9.47076 

298 

90 

A10 

40.00000 

16.17731 

14.64463 

322 

97 

All 

43.00000 

16.39803 

14.14314 

360 

107 

A12 

63.00000 

10.03833 

9.76055 

403 

7 

A13 

133.00000 

17.36060 

12.01885 

382 

32 

A14 

69.00000 

10.67624 

7.85616 

417 

32 

A15 

62.00000 

13.89731 

10.04098 

404 

40 

A16 

59.00000 

13.12808 

10.83174 

474 

8 

A17 

134.00000 

15.68090 

12.23588 

506 

149 

A18 

47.00000 

15.95534 

13.42601 

508 

171 

A19 

39.00000 

16.21332 

14.80555 

475 

140 

A20 

72.00000 

9.27506 

6.72040 

463 

160 

A21 

154.00000 

16.62390 

15.72015 

454 

182 

A22 

167.00000 

15.79108 

17.41017 

441 

203 

A23 

76.00000 

7.53977 

4.93013 

432 

204 

A24 

164.00000 

17.55718 

17.04832 

428 

211 

A25 

87.00000 

-1.42042 

-3.69298 

393 

155 

A26 

75.00000 

6.84302 

5.41346 

377 

216 

A27 

167.00000 

16.23982 

17.41017 

346 

250 

A28 

122.00000 

10.70919 

9.18588 

285 

207 

A29 

63.00000 

12.31237 

9.76055 

262 

192 

A30 

63.00000 

10.28495 

9.76055 

418 

280 

A31 

117.00000 

5.72776 

7.51193 

420 

311 

A32 

70.00000 

9.19157 

7.49383 

440 

302 

A33 

50.00000 

14.35641 

12.84642 

478 

256 

A34 

117.00000 

7.50255 

7.51193 

454 

392 

A35 

126.00000 

8.33095 

10.32574 

492 

395 

A36 

56.00000 

13.44685 

11.55648 

489 

401 

A37 

127.00000 

9.00575 

10.58877 

478 

498 

A38 

89.00000 

-9.69323 

-6.90161 

9 

330 

A39 

119.00000 

6.50417 

8.22067 

98 

291 

A40 

44.00000 

12.83869 

13.96933 

<18 

303 

A41 

64.00000 

9.10144 

9.47076 

170 

341 

A42 

115.00000 

10.04106 

6.74018 

121 

373 

A43 

67.00000 

9.16296 

8.53842 

183 

434 

A44 

112.00000 

-3.26578 

5.43645 

198 

435 

A45 

124.00000 

8.41299 

9.77449 

272 

413 

A46 

56.00000 

12.40393 

11.55648 
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Table  C-2.  Channel  5,  no  preprocessing. 


Line 

Pixel 

133 

50 

A1A 

96.00000 

14.89897 

14.52934 

133 

51 

A1B 

94.00000 

15.27888 

14.90951 

137 

76 

A2A 

97.00000 

15.74141 

14.33283 

137 

77 

A2B 

70.00000 

19.94801 

18.51143 

162 

51 

A3A 

118.00000 

9.78546 

8.64409 

163 

50 

A3B 

113.00000 

11.36008 

10.37303 

57 

194 

A4A 

94.00000 

14.41450 

14.90951 

58 

195 

A4B 

85.00000 

15.79175 

16.44022 

237 

66 

A5A 

80.00000 

18.42231 

17.18686 

237 

67 

A5B 

99.00000 

15.11119 

13.92594 

245 

158 

A6A 

148.00000 

0.48357 

-1.21316 

246 

158 

A6B 

162.00000 

9.86171 

8.08956 

199 

92 

A7 

172.00000 

10.05867 

11.43312 

188 

235 

A8 

164.00000 

7.80173 

8.86916 

259 

119 

A9 

97.00000 

15.38368 

14.33283 

298 

90 

A10 

53.00000 

19.34751 

20.38299 

322 

97 

All 

33.00000 

20.47749 

22.16727 

360 

107 

A12 

75.00000 

19.52915 

17.87437 

403 

7 

A13 

188.00000 

16.25512 

15.01965 

382 

32 

A14 

101.00000 

13.90491 

13.49905 

417 

32 

A15 

82.00000 

16.34279 

16.89586 

404 

40 

A16 

69.00000 

17.22554 

18.63342 

474 

8 

A17 

194.00000 

,  16.62584 

16.05701 

506 

149 

A18 

54.00000 

18.84574 

20.28338 

508 

171 

A19 

49.00000 

20.14812 

20.77038 

475 

140 

A20 

95.00000 

14.24431 

14.72151 

463 

160 

A21 

213.00000 

16.37731 

18.70532 

454 

182 

A22 

190.00000 

12.12395 

15.37938 

441 

203 

A23 

102.00000 

11.26938 

13.27747 

432 

204 

A24 

197.00000 

14.43143 

16.53265 

428 

211 

A25 

122.00000 

4.90524 

6.96057 

393 

155 

A26 

138.00000 

-13.12124 

-9.84341 

377 

216 

A27 

223.00000 

17.18371 

19.83058 

346 

250 

A28 

160.00000 

5.10370 

7.23303 

285 

207 

A29 

96.00000 

14.40050 

14.52934 

262 

192 

A30 

93.00000 

13.03213 

15.09353 

418 

280 

A31 

133.00000 

-1.03939 

-0.74510 

420 

311 

A32 

79.00000 

15.60697 

17.32878 

440 

302 

A33 

77.00000 

16.82623 

17.60586 

478 

256 

A34 

138.00000 

-11.40859 

-9.84341 

454 

392 

A35 

152.00000 

0.36730 

2.58281 

492 

395 

A36 

65.00000 

16.63120 

19.10495 

489 

401 

A37 

160.00000 

4.95219 

7.23303 

478 

498 

A38 

144.00000 

-10.54844 

-8.11226 

9 

330 

A39 

151.00000 

-1.88051 

1.77783 

98 

291 

A40 

25.00000 

21.08889 

22.78962 

118 

303 

A41 

78.00000 

15.47009 

17.46843 

170 

341 

A42 

162.00000 

8.64341 

8.08956 

121 

373 

A43 

85.00000 

13.50184 

16.44022 

183 

434 

A44 

131.00000 

-0.26179 

1.26005 

198 

435 

A45 

150.00000 

-0.80553 

0.89056 

272 

413 

A46 

58.00000 

16.75670 

19.87312 
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Table  C-3.  Channel  1,  preprocessed. 


Line 

133 

50 

A1A 

87.00000 

13.18351 

12.39540 

133 

51 

A1B 

92.00000 

11.99214 

11.17592 

137 

76 

A2A 

90.00000 

13.72329 

11.68443 

137 

77 

A2B 

62.00000 

18.80811 

16.77104 

162 

51 

A3A 

108.00000 

6.70344 

5.45678 

163 

50 

A3B 

102.00000 

9.15027 

8.05970 

57 

194 

A4A 

84.00000 

12.56092 

13.05255 

58 

195 

A4B 

81.00000 

12.96470 

13.66346 

237 

66 

A5A 

77.00000 

16.00595 

14.41642 

237 

67 

A5B 

90.00000 

13.20043 

11.68443 

245 

158 

A6A 

122.00000 

-7.34266 

-9.21409 

246 

158 

A6B 

135.00000 

2.35190 

0.61036 

199 

92 

A7 

165.00000 

11.24455 

12.76449 

188 

235 

A8 

135.00000 

-0.53703 

0.61036 

259 

119 

A9 

91.00000 

12.86638 

11.43389 

298 

90 

A10 

50.00000 

17.92134 

18.28168 

322 

97 

All 

31.00000 

18.97251 

20.23900 

360 

107 

A12 

72.00000 

17.65183 

15.27419 

403 

7 

A13 

185.00000 

17.86539 

16.29877 

382 

32 

A14 

94.00000 

11.54459 

10.63577 

417 

32 

A15 

70.00000 

15.79424 

15.59491 

404 

40 

A16 

63.00000 

15.86881 

16.63244 

474 

8 

A17 

192.00000 

18.07628 

17.25985 

506 

149 

A18 

56.00000 

16.68670 

17.55916 

508 

171 

A19 

51.00000 

18.36644 

18.16535 

475 

140 

A20 

89.00000 

11.64286 

11.92794 

463 

160 

A21 

217.00000 

17.64830 

20.02016 

454 

182 

A22 

185.00000 

12.96514 

16.29877 

441 

203 

A23 

95.00000 

8.28025 

10.35255 

432 

204 

A24 

201.00000 

16.21811 

18.35726 

428 

211 

A25 

112.00000 

1.09238 

3.15332 

393 

155 

A26 

122.00000 

-11.46995 

-9.21409 

377 

216 

A27 

223.00000 

17.88440 

20.56992 

346 

250 

A28 

154.00000 

7.85478 

9.95705 

285 

207 

A29 

91.00000 

11.99037 

11.43389 

262 

192 

A30 

83.00000 

11.25739 

13.26099 

418 

280 

A31 

130.00000 

-6.61097 

-5.56314 

420 

311 

A32 

78.00000 

12.99186 

14.23421 

440 

302 

A33 

73.00000 

14.81421 

15.10928 

478 

256 

A34 

144.00000 

5.54500 

6.25693 

454 

392 

A35 

151.00000 

6.36485 

9.00260 

492 

395 

A3  6 

55.00000 

15.00566 

17.68383 

489 

401 

A37 

158.00000 

8.60537 

11.08563 

478 

498 

A38 

129.00000 

-9.33133 

-7.57882 

9 

330 

A39 

135.00000 

-3.37199 

0.61036 

98 

291 

A40 

23.00000 

19.80061 

20.94721 

118 

303 

A41 

75.00000 

13.40707 

14.76975 

170 

341 

A42 

145.00000 

7.00428 

6.70641 

121 

373 

A43 

'76.00000 

12.10294 

14.59488 

183 

434 

A44 

126.00000 

-22.26839 

-20.87355 

198 

435 

A45 

144.00000 

4.53005 

6.25693 

272 

413 

A46 

53.00000 

15.09171 

17.92793 

C-5 


Table  C-4.  Channel  5,  preprocessed. 


133 

50 

A1A 

96.00000 

14.89897 

14.52934 

133 

51 

A1B 

94.00000 

15.27888 

14.90951 

137 

76 

A2A 

97.00000 

15.74141 

14.33283 

137 

77 

A2B 

70.00000 

19.94801 

18.51143 

162 

51 

A3A 

118.00000 

9.78546 

8.64409 

163 

50 

A3B 

113.00000 

11.36008 

10.37303 

57 

194 

A4A 

94.00000 

14.41450 

14.90951 

58 

195 

A4B 

85.00000 

15.79175 

16.44022 

237 

66 

A5A 

80.00000 

18.42231 

17.18686 

237 

67 

A5B 

99.00000 

15.11119 

13.92594 

245 

158 

A6A 

148.00000 

0.48357 

-1.21316 

246 

158 

A6B 

162.00000 

9.86171 

8.08956 

199 

92 

A7 

172.00000 

10.05867 

11.43312 

188 

235 

A8 

164.00000 

7.80173 

8.86916 

259 

119 

A9 

97.00000 

15.38368 

14.33283 

298 

90 

A10 

53.00000 

19.34751 

20.38299 

322 

97 

All 

33.00000 

20.47749 

22.16727 

360 

107 

A12 

75.00000 

19.52915 

17.87437 

403 

7 

A13 

188.00000 

16.25512 

15.01965 

382 

32 

A14 

101.00000 

13.90491 

13.49905 

417 

32 

A15 

82.00000 

16.34279 

16.89586 

404 

40 

A16 

69.00000 

17.22554 

18.63342 

474 

8 

A17 

194.00000 

16.62584 

16.05701 

506 

149 

A18 

54.00000 

18.84574 

20.28338 

508 

171 

A19 

49.00000 

20.14812 

20.77038 

475 

140 

A20 

95.00000 

14.24431 

14.72151 

463 

160 

A21 

213.00000 

16.37731 

18.70532 

454 

182 

A22 

190.00000 

12.12395 

15.37938 

441 

203 

A23 

102.00000 

11.26938 

13.27747 

432 

204 

A24 

197.00000 

14.43143 

16.53265 

428 

211 

A25 

122.00000 

4.90524 

6.96057 

393 

155 

A26 

138.00000 

-13.12124 

-9.84341 

377 

216 

A27 

223.00000 

17.18371 

19.83058 

346 

250 

A28 

160.00000 

5.10370 

7.23303 

285 

207 

A29 

96.00000 

14.40050 

14.52934 

262 

192 

A30 

93.00000 

13.03213 

15.09353 

418 

280 

A31 

133.00000 

-1.03939 

-0.74510 

420 

311 

A32 

79.00000 

15.60697 

17.32878 

440 

302 

A33 

77.00000 

16.82623 

17.60586 

478 

256 

A34 

138.00000 

-11.40859 

-9.84341 

454 

392 

A35 

152.00000 

0.36730 

2.58281 

492 

395 

A3  6 

65.00000 

16.63120 

19.10495 

489 

401 

A37 

160.00000 

4.95219 

7.23303 

478 

498 

A38 

144.00000 

-10.54844 

-8.11226 

9 

330 

A39 

151.00000 

-1.88051 

1.77783 

98 

291 

A40 

25.00000 

21.08889 

22.78962 

118 

303 

A41 

78.00000 

15.47009 

17.46843 

170 

341 

A42 

162.00000 

8.64341 

8.08956 

121 

373 

A43 

85.00000 

13.50184 

16.44022 

183 

434 

A44 

131.00000 

-0.26179 

1.26005 

198 

435 

A45 

150.00000 

-0.80553 

0.89056 

272 

413 

A46 

58.00000 

16.75670 

19.87312 

C-6 


Table  C-5.  Minimum  noise,  channels  1  and  5,  no  preprocessing. 


Line 

Local  SNR 

133 

50 

A1A 

31.00000 

15.38540 

15.36220 

133 

51 

A1B 

36.00000 

19.46432 

19.96654 

137 

76 

A2A 

35.00000 

19.01599 

19.22065 

137 

77 

A2B 

36.00000 

19.92917 

19.96654 

162 

51 

A3A 

28.00000 

14.42188 

10.63891 

163 

50 

A3B 

28.00000 

14.43850 

10.63891 

57 

194 

A4A 

33.00000 

17.83999 

17.50394 

58 

195 

A4B 

36.00000 

20.49207 

19.96654 

237 

66 

A5A 

35.00000 

20.54815 

19.22065 

237 

67 

A5B 

31.00000 

16.87674 

15.36220 

245 

158 

A6A 

14.00000 

19.90649 

18.14034 

246 

158 

A6B 

13.00000 

20.90228 

18.98038 

199 

92 

A7 

27.00000 

11.34620 

8.24558 

188 

235 

A8 

15.00000 

17.34940 

17.21027 

259 

119 

A9 

32.00000 

18.48232 

16.49891 

298 

90 

A10 

38.00000 

20.81809 

21.28993 

322 

97 

All 

38.00000 

21.70469 

21.28993 

360 

107 

A12 

33.00000 

18.79750 

17.50394 

403 

7 

A13 

28.00000 

15.85151 

10.63891 

382 

32 

A14 

30.00000 

17.24889 

14.05400 

417 

32 

A15 

28.00000 

13.00053 

10.63891 

404 

40 

A16 

33.00000 

17.90347 

17.50394 

474 

8 

A17 

29.00000 

14.26094 

12.51323 

506 

149 

A18 

36.00000 

12.17385 

19.96654 

508 

171 

A19 

35.00000 

12.32998 

19.22065 

475 

140 

A20 

31.00000 

17.8592C 

15.36220 

463 

160 

A21 

33.00000 

19.47829 

17.50394 

454 

182 

A22 

30.00000 

15.31729 

14.05400 

441 

203 

A23 

32.00000 

16.30754 

16.49891 

432 

204 

A24 

36.00000 

19.45976 

19.96654 

428 

211 

A25 

26.00000 

5.50052 

4.93086 

393 

155 

A26 

32.00000 

10.96959 

16.49891 

377 

216 

A27 

29.00000 

13.74139 

12.51323 

346 

250 

A28 

31.00000 

17.98290 

15.36220 

285 

207 

A29 

35.00000 

19.60094 

19.22065 

262 

192 

A30 

30.00000 

13.99266 

14.05400 

418 

280 

A31 

28.00000 

7.30920 

10.63891 

420 

311 

A32 

36.00000 

18.22665 

19.96654 

440 

302 

A33 

37.00000 

18.25457 

20.65342 

478 

256 

A34 

36.00000 

18.75273 

19.96654 

454 

392 

A35 

36.00000 

20.42850 

19.96654 

492 

395 

A3  6 

33.00000 

17.09636 

17.50394 

489 

401 

A37 

33.00000 

17.33883 

17.50394 

478 

498 

A38 

24.00000 

-5.41655 

-18.08220 

9 

330 

A39 

26.00000 

1.53669 

4.93086 

98 

291 

A40 

41.00000 

22.11342 

22.95985 

118 

303 

A41 

37.00000 

18.09341 

20.65342 

170 

341 

A42 

20.00000 

7.67877 

9.97790 

121 

373 

A43 

32.00000 

6.49315 

16.49891 

183 

434 

A44 

29.00000 

11.08483 

12.51323 

198 

435 

A45 

29.00000 

12.37731 

12.51323 

272 

413 

A46 

35.00000 

17.13851 

19.22065 

C-7 


Table  C-6.  Minimum  noise,  channels  5  and  1,  no  preprocessing. 


Line 

Pixel 

Name 

M  . . 

133 

50 

A1A 

26.00000 

15.37291 

15.53756 

133 

51 

A1B 

21.00000 

19.45158 

20.11377 

137 

76 

A2A 

22.00000 

19.01599 

19.37149 

137 

77 

A2B 

21.00000 

19.92917 

20.11377 

162 

51 

A3A 

29.00000 

14.42188 

10.86345 

163 

50 

A3B 

29.00000 

14.43850 

10.86345 

57 

194 

A4A 

24.00000 

18.06302 

17.66436 

58 

195 

A4B 

21.00000 

20.80436 

20.11377 

237 

66 

ASA 

22.00000 

20.91653 

19.37149 

237 

67 

A5B 

26.00000 

17.25636 

15.53756 

245 

158 

A6A 

43.00000 

19.88492 

18.19705 

246 

158 

A6B 

44.00000 

20.88132 

19.04172 

199 

92 

A7 

32.00000 

3.73474 

0.02730 

188 

235 

A8 

43.00000 

18.16886 

18.19705 

259 

119 

A9 

25.00000 

18.48232 

16.66589 

298 

90 

A10 

19.00000 

20.87400 

21.43145 

322 

97 

All 

19.00000 

21.89744 

21.43145 

360 

107 

A12 

24.00000 

19.03107 

17.66436 

403 

7 

A13 

31.00000 

12.51548 

5.26339 

382 

32 

A14 

27.00000 

17.24889 

14.24045 

417 

32 

A15 

29.00000 

12.97928 

10.86345 

404 

40 

A16 

24.00000 

17.95919 

17.66436 

474 

8 

A17 

30.00000 

11.37887 

8.50713 

506 

149 

A18 

21.00000 

12.13800 

20.11377 

508 

171 

A19 

22.00000 

12.33486 

19.37149 

475 

140 

A20 

26.00000 

17.67335 

15.53756 

463 

160 

A21 

26.00000 

17.84381 

15.53756 

454 

182 

A22 

29.00000 

12.18016 

10.86345 

441 

203 

A23 

25.00000 

16.73686 

16.66589 

432 

204 

A24 

23.00000 

18.33480 

18.55980 

428 

211 

A25 

31.00000 

5.96738 

5.26339 

393 

155 

A26 

25.00000 

10.94818 

16.66589 

377 

216 

A27 

30.00000 

11.26104 

8.50713 

346 

250 

A28 

28.00000 

15.75177 

12.71504 

285 

207 

A29 

22.00000 

19.56598 

19.37149 

262 

192 

A30 

27.00000 

14.00073 

14.24045 

418 

280 

A31 

30.00000 

4.86304 

8.50713 

420 

311 

A32 

21.00000 

18.10599 

20.11377 

440 

302 

A33 

20.00000 

17.94036 

20.79757 

478 

256 

A34 

21.00000 

18.40151 

20.11377 

454 

392 

A35 

23.00000 

20.30154 

18.55980 

492 

395 

A36 

24.00000 

18.75406 

17.66436 

489 

401 

A37 

26.00000 

16.66735 

15.53756 

478 

498 

A38 

33.00000 

-4.94467 

-15.22566 

9 

330 

A39 

33.00000 

-6.13650 

-15.22566 

98 

291 

A40 

16.00000 

21.42084 

23.09532 

118 

303 

A41 

20.00000 

17.43640 

20.79757 

170 

341 

A42 

39.00000 

11.42710 

13.63571 

121 

373 

A43 

25.00000 

6.53367 

16.66589 

183 

434 

A44 

28.00000 

11.07893 

12.71504 

198 

435 

A45 

30.00000 

8.10274 

8.50713 

272 

413 

A46 

22.00000 

17.04173 

19.37149 

08 


Table  C-7.  Minimum  noise,  channels  1  and  5,  preprocessed. 


IrfMB 

m  ii  mil 

133 

50 

A1A 

50.00000 

16.02261 

17.45886 

133 

51 

A1B 

61.00000 

20.32410 

22.11472 

137 

76 

A2A 

53.00000 

17.98388 

18.99471 

137 

77 

A2B 

54.00000 

18.65380 

19.45172 

162 

51 

A3A 

46.00000 

16.14147 

14.86821 

163 

50 

A3B 

45.00000 

15.36251 

14.07879 

57 

194 

A4A 

50.00000 

17.03802 

17.45886 

58 

195 

A4B 

59.00000 

21.14482 

21.43340 

237 

66 

A5A 

61.00000 

21.17553 

22.11472 

237 

67 

A5B 

51.00000 

17.15791 

18.00155 

245 

158 

A6A 

18.00000 

18.52890 

17.99056 

246 

158 

A6B 

14.00000 

20.52005 

19.87703 

199 

92 

A7 

43.00000 

12.18818 

12.24541 

188 

235 

A8 

12.00000 

18.76796 

20.68600 

259 

119 

A9 

54.00000 

19.91112 

19.45172 

298 

90 

A10 

65.00000 

20.31334 

23.33534 

322 

97 

All 

69.00000 

21.86564 

24.40539 

360 

107 

A12 

61.0^00 

21.22831 

22.11472 

403 

7 

A13 

46.0*,  00 

14.78591 

14.86821 

382 

32 

A14 

53.00000 

18.19684 

18.99471 

417 

32 

A15 

47.00000 

12.96003 

15.59183 

404 

40 

A16 

59.00000 

18.03839 

21.43340 

474 

8 

A17 

48.00000 

14.54017 

16.25976 

506 

149 

A18 

72.00000 

18.69081 

25.12942 

508 

171 

A19 

73.00000 

19.77211 

25.35794 

475 

140 

A20 

55.00000 

19.66318 

19.88588 

463 

160 

A21 

54.00000 

18.93725 

19.45172 

454 

182 

A22 

45.00000 

12.47961 

14.07879 

441 

203 

A23 

54.00000 

16.68949 

19.45172 

432 

204 

A24 

56.00000 

17.07970 

20.29936 

428 

211 

A25 

45.00000 

10.78576 

14.07879 

393 

155 

A26 

34.00000 

-17.51077 

-12.55745 

377 

216 

A27 

46.00000 

15.04015 

14.86821 

346 

250 

A28 

48.00000 

17.13220 

16.25976 

285 

207 

A29 

57.00000 

18.86016 

20.69405 

262 

192 

A30 

50.00000 

15.70276 

17.45886 

418 

280 

A31 

55.00000 

14.27230 

19.88588 

420 

311 

A32 

64.00000 

19.32911 

23.04589 

440 

302 

A33 

60.00000 

17.29492 

21.78074 

478 

256 

A34 

67.00000 

20.04259 

23.88683 

454 

392 

A35 

54.00000 

19.06127 

19.45172 

492 

395 

A36 

53.00000 

17.69260 

18.99471 

489 

401 

A37 

53.00000 

17.51050 

18.99471 

478 

498 

A38 

36.00000 

-7.06374 

-2.77157 

9 

330 

A39 

33.00000 

-3.62896 

-2.89245 

98 

291 

A40 

71.00000 

22.43418 

24.89472 

118 

303 

A41 

61.00000 

18.32983 

22.11472 

170 

341 

A42 

29.00000 

4.14606 

8.43713 

121 

373 

A43 

51.00000 

9.79541 

18.00155 

183 

434 

A44 

52.00000 

16.09825 

18.51232 

198 

435 

A45 

47.00000 

13.81605 

15.59183 

272 

413 

A46 

62.00000 

18.78494 

22.43633 

C-9 


Table  C-8.  Minimum  noise,  channels  5  and  1;  preprocessed. 


Line 

Pixel 

Name 

BIkH 

Global  SNR 

133 

50 

A1A 

29.00000 

16.65678 

18.38717 

133 

51 

A1B 

18.00000 

20.54447 

22.62886 

137 

76 

A2A 

26.00000 

18.38977 

19.76359 

137 

77 

A2B 

24.00000 

19.40650 

20.57326 

162 

51 

A3A 

34.00000 

16.57766 

15.45886 

163 

50 

A3B 

35.00000 

15.86384 

14.73286 

57 

194 

A4A 

30.00000 

17.49289 

17.87522 

58 

195 

A4B 

21.00000 

21.42741 

21.66173 

237 

66 

A5A 

18.00000 

21.59146 

22.62886 

237 

67 

A5B 

29.00000 

17.42889 

18.38717 

245 

158 

A6A 

63.00000 

18.54029 

17=90567 

246 

158 

A6B 

67.00000 

20.52421 

19.78810 

199 

92 

A7 

39.00000 

11.26801 

11.00844 

188 

235 

A8 

70.00000 

19.06575 

20.97282 

259 

119 

A9 

25.00000 

20.32002 

20.17785 

298 

90 

A10 

13.00000 

20.65070 

24.03417 

322 

97 

All 

8.00000 

22.38051 

25.24347 

360 

107 

A12 

18.00000 

21.39224 

22.62886 

403 

7 

A13 

36.00000 

13.70128 

13.94061 

382 

32 

A14 

27.00000 

18.29110 

19.32858 

417 

32 

A15 

31.00000 

14.37622 

17.33119 

404 

40 

A16 

20.00000 

18.34177 

21.99621 

474 

8 

A17 

35.00000 

12.86880 

14.73286 

506 

149 

A18 

6.00000 

19.11337 

25.68367 

508 

171 

A19 

5.00000 

20.10514 

25.89568 

475 

140 

A20 

24.00000 

20.14408 

20.57326 

463 

160 

A21 

30.00000 

17.39120 

17.87522 

454 

182 

All 

38.00000 

10.66431 

12.09955 

441 

203 

A23 

26.00000 

17.12608 

19.76359 

432 

204 

A24 

27.00000 

16.22676 

19.32858 

428 

211 

A25 

35.00000 

11.54768 

14.73286 

393 

155 

A26 

46.00000 

-18.14977 

-12.99738 

377 

216 

A27 

38.00000 

12.32922 

12.09955 

346 

250 

A28 

34.00000 

16.22068 

15.45886 

285 

207 

A29 

23.00000 

18.73442 

20.95144 

262 

192 

A30 

29.00000 

16.53267 

18.38717 

418 

280 

A31 

25.00000 

14.57386 

20.17785 

420 

311 

A32 

14.00000 

19.79203 

23.77071 

440 

302 

A33 

19.00000 

17.73445 

22.31829 

478 

256 

A34 

14.00000 

19.86955 

23.77071 

454 

392 

A35 

28.00000 

19.19569 

18.87062 

492 

395 

A36 

25.00000 

19.18844 

20.17785 

489 

401 

A37 

29.00000 

17.26612 

18.38717 

478 

498 

A38 

45.00000 

-7.98254 

-3.10651 

9 

330 

A39 

48.00000 

-2.61790 

-2.77151 

98 

291 

A40 

6.00000 

22.83877 

25.68367 

118 

303 

A41 

18.00000 

18.44745 

22.62886 

170 

341 

A42 

53.00000 

5.64155 

9.83761 

121 

373 

A43 

28.00000 

10.53008 

18.87062 

183 

434 

A44 

29.00000 

15.94966 

18.38717 

198 

435 

A45 

34.00000 

13.67963 

15.45886 

272 

413 

A46 

16.00000 

19.41018 

23.21853 
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Table  C-9.  Simple  difference,  channels  1  and  5,  no  preprocessing. 


Line 

M  II  i  III  II 

133 

50 

A1A 

32.00000 

15.37109 

15.54239 

133 

51 

A1B 

37.00000 

19.45013 

20.12693 

137 

76 

A2A 

36.00000 

19.01599 

19.38358 

137 

77 

A2B 

37.00000 

19.92917 

20.12693 

162 

51 

A3A 

29.00000 

14.42188 

10.85377 

163 

50 

A3B 

29.00000 

14.43850 

10.85377 

57 

194 

A4A 

34.00000 

18.02621 

17.67362 

58 

195 

A4B 

37.00000 

20.73136 

20.12693 

237 

66 

A5A 

36.00000 

20.93970 

19.38358 

237 

67 

A5B 

32.00000 

17.23372 

15.54239 

245 

158 

A6A 

14.00000 

20.61410 

19.08043 

246 

158 

A6B 

14.00000 

20.78595 

19.08043 

199 

92 

A7 

27.00000 

8.51883 

5.22204 

188 

235 

A8 

14.00000 

18.94280 

19.08043 

259 

119 

A9 

33.00000 

18.48232 

16.67321 

298 

90 

A10 

39.00000 

20.87125 

21.44631 

322 

97 

All 

39.00000 

21.84081 

21.44631 

360 

107 

A12 

34.00000 

19.06333 

17.67362 

403 

7 

A13 

28.00000 

14.38342 

8.48656 

382 

32 

A14 

31.00000 

17.30161 

14.24201 

417 

32 

A15 

29.00000 

12.96906 

10.85377 

404 

40 

A16 

34.00000 

17.92350 

17.67362 

474 

8 

A17 

29.00000 

12.86577 

10.85377 

506 

149 

A18 

37.00000 

12.10310 

20.12693 

508 

171 

A19 

36.00000 

12.24652 

19.38358 

475 

140 

A20 

32.00000 

17.78824 

15.54239 

463 

160 

A21 

33.00000 

18.81455 

16.67321 

454 

182 

All 

29.00000 

12.33164 

10.85377 

441 

203 

A23 

33.00000 

16.56220 

16.67321 

432 

204 

A24 

36.00000 

18.94770 

19.38358 

428 

211 

A25 

27.00000 

5.77204 

5.22204 

393 

155 

A26 

33.00000 

10.95981 

16.67321 

377 

216 

All 

29.00000 

13.17847 

10.85377 

346 

250 

A28 

31.00000 

17.19125 

14.24201 

285 

207 

A29 

36.00000 

19.57293 

10.38358 

262 

192 

A30 

31.00000 

13.99416 

14.24201 

418 

280 

A31 

28.00000 

4.90595 

8.48656 

420 

311 

All 

37.00000 

18.13732 

20.12693 

440 

302 

A33 

38.00000 

18.07336 

20.81165 

478 

256 

A34 

37.00000 

18.55390 

20.12693 

454 

392 

A35 

36.00000 

20.76050 

19.38358 

492 

395 

A  36 

34.00000 

18.42007 

17.67362 

489 

401 

All 

33.00000 

17.45537 

16.67321 

478 

498 

All 

25.00000 

-5.06116 

-15.93137 

9 

330 

A39 

26.00000 

-7.69420 

-0.06946 

98 

291 

A40 

42.00000 

21.91499 

23.11197 

118 

303 

A41 

38.00000 

17.77667 

20.81165 

170 

341 

A42 

20.00000 

9.74230 

12.04488 

121 

373 

A43 

33.00000 

6.52587 

16.67321 

183 

434 

A44 

29.00000 

9.60797 

10.85377 

198 

435 

A45 

29.00000 

10.69957 

10.85377 

272 

413 

A46 

36.00000 

17.04091 

19.38358 
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Table  C-10.  Simple  difference,  channels  1  and  5,  preprocessed. 


Line 

Pixel 

Local  SNR 

133 

50 

A1A 

51.00000 

16.22560 

17.78472 

133 

51 

A1B 

62.00000 

20.34799 

22.27235 

137 

76 

A2A 

54.00000 

18.05914 

19.25525 

137 

77 

A2B 

56.00000 

19.12528 

20.11330 

162 

51 

A3A 

47.00000 

16.62221 

15.33251 

163 

50 

A3B 

45.00000 

15.06612 

13.78596 

57 

194 

A4A 

51.00000 

17.40089 

17.78472 

58 

195 

A4B 

60.00000 

21.38097 

21.61055 

237 

66 

A5A 

62.00000 

21.31522 

22.27235 

237 

67 

A5B 

52.00000 

17.43573 

18.30306 

245 

158 

A6A 

18.00000 

18.63743 

18.03507 

246 

158 

A6B 

14.00000 

20.60875 

19.89636 

199 

92 

A7 

43.00000 

12.04107 

11.90307 

188 

235 

A8 

12.00000 

18.79436 

20.69607 

259 

119 

A9 

55.00000 

20.02104 

19.69487 

298 

90 

A10 

67.00000 

20.52220 

23.73501 

322 

97 

All 

72.00000 

22.21943 

24.98650 

360 

107 

A12 

63.00000 

21.50818 

22.58525 

403 

7 

A13 

46.00000 

14.45364 

14.59361 

382 

32 

A14 

54.00000 

18.38504 

19.25525 

417 

32 

A15 

49.00000 

13.84769 

16.64487 

404 

40 

A16 

61.00000 

18.42846 

21.94775 

474 

8 

A17 

47.00000 

13.54165 

15.33251 

506 

149 

A18 

74.00000 

18.89304 

25.44062 

508 

171 

A19 

75.00000 

19.92278 

25.65907 

475 

140 

A20 

56.00000 

19.75880 

20.11330 

463 

160 

A21 

53.00000 

18.32460 

18.79220 

454 

182 

A22 

44.00000 

11.43797 

12.89544 

441 

203 

A23 

55.00000 

17.02972 

19.69487 

432 

204 

A24 

56.00000 

16.97133 

20.11330 

428 

211 

A25 

46.00000 

11.35706 

14.59361 

393 

155 

A26 

35.00000 

-23.28901 

-18.06498 

377 

216 

A27 

45.00000 

14.07800 

13.78596 

346 

250 

A28 

48.00000 

16.87482 

16.01345 

285 

207 

A29 

58.00000 

18.84143 

20.89415 

262 

192 

A30 

51.00000 

15.99623 

17.78472 

418 

280 

A31 

56.00000 

14.47204 

20.11330 

420 

311 

A32 

66.0C000 

19.64910 

23.46152 

440 

302 

A33 

62.00000 

17.69850 

22.27235 

478 

256 

A34 

68.00000 

20.10480 

24.00015 

454 

392 

A35 

54.00000 

19.24547 

19.25525 

492 

395 

A3  6 

55.00000 

18.56292 

19.69487 

489 

401 

A37 

53.00000 

17.47102 

18.79220 

478 

498 

A38 

36.00000 

-8.99934 

-4.41606 

9 

330 

A39 

33.00000 

-1.87627 

-1.63907 

98 

291 

A40 

74.00000 

22.77999 

25.44062 

118 

303 

A41 

62.00000 

18.26562 

22.27235 

170 

341 

A42 

29.00000 

4.44822 

8.73597 

121 

373 

A43 

53.00000 

10.47408 

18.79220 

183 

434 

A44 

52.00000 

15.91463 

18.30306 

198 

435 

A45 

47.00000 

13.53255 

15.33251 

272 

413 

A46 

63.00000 

18.88253 

22.58525 
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APPENDIX  D:  DETAILS  OF  NUMERIC  RESULTS, 
HUNTSVILLE  DATA  SET 


The  TIMS  Huntsville  data  set  was  processed  to  confirm  results  achieved  by  using  the  TIMS 
Adelaide  data  set.  The  TIMS  Adelaide  data  set  was  processed  to  compare  the  results  to  the  spatial 
difference  model  developed  by  L.  E.  Hoff. 

The  first  step  in  this  processing  was  to  determine  which  two  of  the  six  TIMS  bands  met  the 
criterion  of  high  correlation  between  the  bands.  A  routine  was  run  on  the  raw  data  that  produces  a 
correlation  coefficient  matrix,  as  shown  in  figure  D-l.  It  was  also  desirable  for  this  processing  that  the 
bands  be  widely  separated,  which  tends  to  preserve  target  coloring.  Bands  1  and  4,  with  a  correlation 
coefficient  of  0.9015,  were  chosen. 


1.000*10° 

9.750*10-1 

9.376*101 

9.015*10*1 

8.300*10-1 

7.647*10-1 

\ 

9.750*10-1 

1.000*10° 

9.374*10-1 

8.827*10-1 

7.938*10*1 

7.315*10-1 

9.376*10-1 

9.374*101 

1.000*10° 

8.034*10-1 

7.314*10-1 

5.921*10*1 

9.015*10-1 

8.827*10-1 

8.034*10-1 

1.000*10° 

9.625*10-1 

8.958*10-1 

8.300*10-1 

7.938*10-1 

7.314*10-1 

9.625*10-1 

1.000*10° 

9.084*10-1 

/ 

7.647*10-1 

7.315*10-1 

5.921*10-1 

8.958*10-1 

9.084*10-1 

1.000*10° 

Figure  D-l.  Correlation  coefficient  matrix. 


The  second  step  in  the  processing  was  to  window  the  area  of  interest.  An  area  that  was 
256  X  256  in  size  was  chosen.  From  this  area,  22  target  pixels  were  chosen  for  this  processing  run. 
Figure  D-2  identifies  the  target  pixels  by  line  and  column. 


Line 

Pixel 

Pixel  Name 

Line 

Pixel 

Pixel  Name 

160 

114 

TAl 

216 

189 

R2 

159 

113 

TA2 

216 

190 

R3 

159 

114 

TA3 

216 

191 

R4 

104 

181 

TCI 

216 

192 

R5 

104 

182 

TC2 

216 

193 

R6 

103 

180 

TC3 

217 

188 

R7 

103 

181 

TC4 

217 

189 

R8 

56 

148 

SI 

217 

190 

R9 

56 

149 

S2 

217 

191 

R10 

56 

150 

S3 

217 

192 

Rll 

216 

188 

R1 

217 

193 

R12 

Figure  D-2.  Huntsville  target  pixels. 


Three  different  preprocessing  methods  were  used  on  this  data  set.  The  first  case  did  a  global 
normalization  on  the  image,  which  took  the  mean  and  variance  of  the  entire  256  X  256  image  to 
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compute  pixel  values.  The  second  case  did  a  line-by-line  normalization  on  the  data,  which  entailed 
taking  the  mean  and  variance  of  each  line  to  determine  the  pixel  values  and  then  doing  a  global 
normalization  on  the  data.  The  third  case  did  a  line-by-line  normalization  on  the  data,  followed  by  a 
high-pass  5x5  filter.  As  a  last  step,  the  data  set  was  globally  normalized.  The  final  output  of  the 
preprocessing  phase  of  this  image  processing  consisted  of  six  files,  three  each  of  band  1  and  band  4, 
globally  normalized,  line-by-line  normalized,  and  filtered. 

For  each  case,  and  at  each  target  pixel,  the  band  with  the  largest  signal-to-noise  (SNR)  was 
determined  and  the  data  set  was  entered  into  a  database  (see  tables  D-l,  D-2,  and  D-3).  The  pixel 
values  that  determined  the  color  ratio  for  comparison  purposes  were  also  entered  at  this  time. 

Upon  completion  of  the  preprocessing  steps,  three  weighted-difference  algorithms  were  applied  to 
each  set  (band  1  and  band  4)  of  data.  Each  of  these  algorithms  used  a  different  weighting  factor;  one 
maximized  the  output  SNR,  Max_SNR;  one  minimized  the  noise,  Min_Noise;  and  one  used  a  weight¬ 
ing  factor  of  1,  Sim_Diff.  The  Min_Noise  algorithm  was  applied  twice  to  each  data  set.  One  run 
applied  band  1  to  band  4;  the  second  run  applied  band  4  to  band  1.  The  SNR  for  each  target  pixel 
from  each  output  image  was  then  measured  and  entered  into  the  database  with  the  input  band’s  SNRs 
(tables  D-l  through  D-3).  The  dB  gain  was  then  computed  from  the  input  and  output  SNR  informa¬ 
tion  and  entered  into  the  database.  The  data  provided  the  necessary  information  for  comparing  the 
spatial-difference  model  to  the  Huntsville  data. 
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Table  D-l.  Huntsville  normalized  (entire  image  only)  channels  1  and  4  target  data. 
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Table  D-2.  Huntsville  normalized  (line  by  line,  then  entire  image)  channels  1  and  4  target  data. 
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Table  D-3.  Huntsville  normalized  (line  by  line/high-pass  filter/then  renormalized)  channels  1  and  4  target  data. 
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